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A novel Bloch-waves based one-step theory of photoemission is de- 
veloped within the augmented plane wave formalism. Implications of 
multi-Bloch-wave structure of photoelectron final states for band map- 
ping are established. Interference between Bloch components of initial 
and final states leads to prominent spectral features with characteristic 
frequency dispersion experimentally observed in VSe2 and TiTe2. In- 
terference effects together with a non-free-electron nature of final states 
strongly limit the applicability of the common direct transitions band 
mapping approach, making the tool of one-step analysis indispensable. 

Angle resolved photoemission (ARPES) is the most powerful tool to 
study electronic structure of crystals [1] . The basis for the interpretation of 
experimental results is provided by the one-step theory of photoemission [2], 
in which the photoelectron initial states are one-particle eigenstates of the 
semi-infinite crystal, and the final state is the time reversed LEED state |3>). 
The interaction with the many-electron system of the crystal gives rise to 
inelastic scattering and leads to a finite mean free path of the photoelectron. 
In spite of the lacking translational invariance in the surface perpendicular 
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direction the concept of 3-dimensional band structure has proved very fruit- 
ful for interpretation of ARPES spectra [3]. In the band structure picture 
the finite photoelectron mean free path is seen as a crystal momentum un- 
certainty of the final state Bloch waves, whose immediate effect is to broaden 
the spectral structures due to direct transitions. For a reasonably small mo- 
mentum broadening it is often assumed that the bulk photoemission can be 
understood in terms of the E(k±) dispersion of occupied and unoccupied 
bands. 

The one-step model, however, implies a more complicated picture: fol- 
lowing Slater [4] and Pendry [5] the outgoing electron is described by a 
pure state represented by a damped wave function: the inelastic scattering 
is taken into account by adding an imaginary part — iV\ to the potential in 
the crystal half-space, so that |3>) is an eigenfunction of a non-Hermitean 
Hamiltonian corresponding to a real eigenvalue Ef - the photoelectron final 
energy. The imaginary term causes the LEED function to decay in space 
and provides the model with the desired surface sensitivity. The final state, 
thus, retains its fully coherent nature, which governs the interference of the 
transition matrix elements between the Bloch constituents of the final and 
the initial state. 

In the simplest case the final state comprises only one dominant Bloch 
wave, as in VSe2 in the energy range shown in Fig. 1, but there may be two 
or more Bloch constituents, as in TiTe2 in Fig. 1. The photoelectron initial 
state is a standing wave, and in the depth of the crystal it is a superposition 
of the Bloch wave |k + ) incident from the interior of the crystal on the surface 
and a number of reflected waves traveling in the opposite direction. In 
the simplest case there is only one reflected wave |k~), see Fig. 1. The 
dipole matrix elements (<I>|p|k + } and ($|p|k _ ) may cancel or enhance each 
other depending on their phases. Let us consider a single-Bloch-wave final 
state. When the absorbing potential is small the phases of the matrix 
elements are unimportant because only the matrix element of the momentum 
conserving transition is large (see the transition at k in the left panel of 
Fig. 1). For realistic values of Vj, which are of the order of 1 eV, the matrix 
elements become comparable, and the interference of the waves becomes 
noticeable. For brevity, we shall refer to this as indirect interference; its 
importance increases with inelastic scattering strength. If the final state 
is composed of several Bloch states strong interference may occur with a 
negligible momentum broadening, as in TiTe2 in Fig. 1. This we shall call 
direct interference. Although a composite structure of final states has been 
noticed earlier [6] , the effect of interference and its implications for the band 
mapping have not been realized. 
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In the present paper we elucidate the influence of the Bloch wave inter- 
ference on the observed spectra with the aim to understand the connection 
between a photoemission spectrum and the underlying band structure in 
more detail. In view of the phenomenological nature of the final state the 
question arises whether the interference effects are correctly treated in the 
one-step theory and whether the multi-Bloch-wave structure of the final 
state may provide additional information about electron states. To clarify 
these issues, we have performed a combined theoretical and experimental 
study of normal emission from VSe2 and TiTe2 layered crystals. 

Our theory is based on the band structure approach both to final and 
to initial states [8]. The LEED state is a solution of a scattering problem 
for a plane wave incident from vacuum. In the bulk of the crystal the wave 
function is given in terms of the complex band structure, and it is continued 
into the surface region by solving a Cauchy problem. A detailed description 
of the methodology for an all-electron potential of general shape within the 
augmented plane wave (APW) formalism has been presented in Ref. [9]. A 
novel aspect of the present method is that the same scattering technique is 
used for initial states, only the incident wave is now the Bloch wave |k + ) 
[10]. 

The APW based Bloch wave scattering theory offers an accurate and 
efficient computational scheme, whose characteristic features are: (i) all 
multiple scattering in the crystal and at the surface is taken into account 
within a full-potential technique; (ii) connection between the spectra and the 
complex band structure of the crystal is straightforward and transparent; 
(iii) no electron absorption is introduced in calculating initial states. The 
treatment of initial states is the main difference between the present the- 
ory and more traditional KKR-based scattering theories: the layer doubling 
method of Ref. [5] essentially relies on the imaginary term in the Hamilto- 
nian, and in real-space techniques [11] the Bloch character of the electron 
states is obscured. 

The ARPES experiment was performed at the SuperACO synchrotron in 
LURE, France. The spectra were measured at normal emission with photon 
energies between 11.5 and 35 eV. The synchrotron radiation polarization 
vector was set at an angle of 45° to the surface normal in the MYM' azimuth. 
The combined monochromator and analyzer energy resolution varied from 
23 to 130 meV with increase of the photon energies through the experimental 
range. 

Experimental and theoretical spectra for VSe2 are compared in Fig. 2. 
Our theory correctly reproduces all the main structures and their disper- 
sion with photon energy. Similar level of agreement is obtained for TiTe2- 



3 



Final states were validated by comparing the energy dependence of the cur- 
rent carried by the LEED state with the electron transmission measurements 
[12]. A detailed comparison of the peak dispersion in theory and experiment 
for VSe2 and TiTe2 is presented in Fig. 3. Well visible are the discrepan- 
cies due to simplified treatment of quasiparticles within the local density 
approximation (LDA): the calculated upper Se Ap z band is shifted upwards 
by 0.3 eV, and the Se 4p xy band downwards by 0.5 eV [13]. Discrepancies 
of the same character and size are observed for the Te 5p states in TiTe2- In 
both crystals, below Ef ~ 16 eV the measured final energy location of the 
peaks is well reproduced, and above 16 eV the measured dispersion curves 
are stretched towards higher energies. It is interesting that although the 
unoccupied bands of VSe2 are shifted by 3 to 4 eV upwards relative to the 
bands of TiTe2 the self-energy effects are seen to start growing at roughly 
the same energy. 

The unoccupied bands look rather complicated, but the calculation points 
out the few waves important for the electron's escape to vacuum. Their 
structure is relatively simple and transparent (in Fig. 3 the conducting bands 
are shown by thick lines in the final-states graphs). We note, however, that 
a nearly free electron model is inapplicable here simply because over wide 
energy regions there are two waves strongly contributing to the LEED func- 
tion. The Bloch-wave scattering theory enables us to explicitly separate 
out the effect of k± dispersion of the initial and final Bloch waves and to 
reveal the more subtle factors due to the structure of the wave functions. 
It is instructive to plot the direct transitions lines - a procedure frequently 
used in comparing experimental spectra with a calculated band structure: 
the dashed lines in Fig. 3 show the points with coordinates [Ef(k±), E[(k±_)], 
with k± running over the ID Brillouin zone. The calculated peak dispersion 
is seen to often deviate from the dashed lines. A vivid example of such an 
intrinsic shift [14] is the downward dispersion of the lower p z bands between 
20 and 21 eV in VSe2 and between 16 and 17 eV in TiTe2- The analysis 
shows that it is not just an overlap of two transitions but their indirect 
interference that determines the dispersion of the peak. 

A clear manifestation of the two-branch structure of the LEED states is 
a curious dispersion of the lower Te 5p z peak between Ef = 11 and 15 eV: it 
disperses downwards but does not reach the bottom of the band and vanishes 
at ~ 13 eV to appear again at the valence band maximum less than 3 eV 
higher in energy. This happens because the leading final state wave hops 
from branch a to branch (3 (see Fig. 1) at Ef « 12 eV. In VSe2 at low energies 
only one Bloch wave dominates (see Fig. 1), which allows us to more clearly 
see the effect: in VSe2 the peak disperses almost to the band minimum, and 
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the maximum is reached some 6 eV higher in energy. Although in TiTe2 the 
direct interference influences the shape of the spectra, it is not essential for 
understanding this behavior: the rapid change from minimum to maximum 
is caused by the k± gap between the a and f3 branches. 

A constructive direct interference plays an important role in formation 
of the structures x in VSe2, around Ef = 24.5 eV at E\ rj —1.2 eV, and in 
TiTe2, around E{ = 21.5 eV at E\ ~ —0.5 eV, see Fig. 3. In both materials 
the experimentally observed intensity enhancement correlates well with the 
crossing points of the two direct transition lines. In VSe2 the structure 
appears well below the valence band maximum, and the wave vectors of 
both final state constituents are well distant from the V point. In TiTe2 the 
dispersion is less pronounced, which is consistent with the differences in the 
band structures (see right panels of Fig. 4): in TiTe2 the crossing occurs 
much closer to the V point than in VSe2, and at the same time the initial 
state dispersion is weaker. 

The interference origin of the features is illustrated by Fig. 4, where 
we compare VSe2 spectra for fiuj = 25 and 26.5 eV calculated within the 
one-step theory with the spectra resulting from incoherent summation of 
the squared moduli of the elements (<I>|p|k + ) and ($|p|k _ ). For a small 
momentum broadening, V\ = 0.25 eV, the deviation of a coherent curve 
from its incoherent counterpart reveals the effect of direct interference. The 
indirect interference leads to differences that are seen in V\ = 1 eV curves 
and not in Vj = 0.25 eV ones. Our important finding is that the structures x 
are due to a strong final-state energy dependence of the direct interference: 
in Fig. 4 the intensity enhancement Dc is large for fiu) = 25 eV and negligible 
for hu = 26.5 eV. 

To illustrate the importance of the interference structures for unambigu- 
ous band mapping, we show in the right panels of Fig. 4 constant initial 
state spectra for the energies E\ at which the strongest intensity enhance- 
ment is observed. For VSe2 the calculated and measured initial energies are 
E[ = —1 and —1.2 eV, respectively (see Fig. 2), and for TiTe2 it is —0.2 and 
—0.5 eV. The final state energy of the direct interference transition is seen to 
be reliably determined in experiment (as expected, the experimental struc- 
tures are shifted to higher Ef energies due to self-energy effects.) Because 
the x structure corresponds to a symmetric k± location of the two final state 
branches the k± assignment of the initial state is stable to the self-energy 
shift of the final states - in contrast to the case of a single-Bloch-wave final 
state. 

To summarize, a Bloch-waves based theory of photoemission is success- 
fully applied to trace back the experimentally observed peak dispersion in 
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VSe2 and TiTe2 to the band structure of initial and final states. Owing 
to Bloch wave interference, information on just dispersion of initial and 
conducting final states is often insufficient to understand the experimental 
results; naive band mapping may be misleading if not supported by accurate 
calculations of intensities. At the same time, our results indicate a domi- 
nant role of direct transitions, which enables reliable band mapping in the 
absence of interference. 

For the first time spectral features are identified that arise from the multi- 
Bloch-wave structure of the final states. Structures of this novel type are 
observed as turning points in the peak dispersion diagrams that correspond 
neither to minima nor to maxima of the initial state bands. Owing to the 
dominant role of direct transitions and the interference enhancement of the 
emission intensity, such points provide information about the wave vectors 
of the final state constituent and allow a more detailed band mapping. 

The deviation of the calculated peak dispersion from the direct transition 
lines is the result of the Bloch wave interference in the initial states. This 
effect is intimately connected to the optical-potential approach to inelastic 
scattering in the one-step theory. The good agreement with experiment in 
cases of strong interference, thus, provides the most direct proof so far in 
favor of the damped waves treatment of inelastic scattering. 

The work was supported by Deutsche Forschungsgemeinschaft via Forscher- 
gruppe FOR 353 and grant CL 124/5-2, and the experiments at LURE by 
the EC within the Access to Research Infrastructure program. 
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Figure 1: Se 4p z and Te 5p z bands of VSe2 and TiTe2 and fragments of 
unoccupied band structure. Full curves display real band structure and 
dashed curves real lines of complex band structure with V\ = 0. Thick 
lines show conducting constituents of the final states calculated with optical 
potential Vj = 1 eV. The line thickness is proportional to the current carried 
(absorbed) by this partial wave [7]. In TiTe2 two conducting branches are 
seen: a and (5. Straight lines illustrate cases of indirect (left panel) and 
direct (right panel) interference. Arrows in the right panel indicate the 
final state k± gap responsible for the dispersion turning point in TiTe2 at 
E{ = 12.5 eV in Fig.3. 
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Figure 2: Normal emission energy distribution curves of VSe2 parametrized 
by photon energy: experiment (left) and theory (right). Spectra are nor- 
malized to the same integral intensity. 
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Figure 3: (color online) Peak dispersion diagrams for normal emission from 
VSe2 and TiTe2. Calculated band structure for initial states is shown in left 
panels and for final states in upper panels with thin lines. In the final states 
panels thick lines show the branches of complex band structure (CBS) that 
most strongly contribute to the final state (thickness is proportional to the 
current absorbed by the Bloch component of the LEED state). In the disper- 
sion diagram dashed lines show the location of direct transitions to the con- 
ducting branches of CBS (lines due to nondispersive p xy and V 3d states are 
not shown). The peak locations in the theoretical EDCs is shown by white 
lines with black fringes (thickness is proportional to peak intensity). They 
are superimposed on the experimental second derivative map —d 2 I{E\)/dEf 
shown in logarithmic colorscale (areas with —d 2 I{E\)/dE? < are clipped). 
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Figure 4: Interference effects in VSe2- Left: Full lines are energy distribution 
curves calculated within the one-step theory, and dashed lines the result 
of incoherent summation of the emission from the |k~) and |k + ) states. 
Structures affected by indirect (I) and direct (D) interference are marked 
by symbols Ic,d and Dc, with subscripts referring to constructive (C) and 
destructive (D) interference. Direct interference is recognized by its influence 
on the spectra with = 0.25 eV. Right: Initial and final states whose direct 
interference causes the feature x in the dispersion diagrams in Fig. 3. Lower 
panels show calculated (lines) and measured (circles) constant initial state 
spectra for the energies E\ at which the strongest intensity enhancement is 
observed. 
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